Alternative genome sequencing approaches of SARS-CoV-2 using Ion AmpliSeq Technology.

A thorough understanding of SARS-CoV-2 genetic features is compulsory to track the ongoing pandemic across multiple geographical locations of the world. Thermo Fisher Scientific USA has developed the Ion AmpliSeq SARS-CoV-2 Research Panel for the targeted sequencing of SARS-CoV-2 complete genome with high coverage and lower error rate. In this study an alternative approach of complete genome sequencing has been validated using different commercial sequencing kits to sequence the SARS-CoV-2. Amplification of cDNA with the SARS-CoV-2 primer pool was performed separately using two different master mixes: 2X environmental master mix (EM) and Platinum™ PCR SuperMix High Fidelity master mix (PM) instead of 5X Ion AmpliSeq™ HiFi Mix whereas NEBNext® Fast DNA Library Prep Set for Ion Torrent™ kit was used as an alternative to Ion AmpliSeq Library Kit Plus for other reagents. This study demonstrated a successful procedure to sequence the SARS-CoV-2 whole genome with average ∼2351 depth and 98.1% of total the reads aligned against the reference sequence (SARS-CoV-2, isolate Wuhan-Hu-1, complete genome). Although genome coverage varied, complete genomes were retrieved for both reagent sets with a reduced cost. This study proposed an alternative approach of high throughput sequencing using Ion torrent technology for the sequencing of SARS-CoV-2 in developing countries where sequencing facilities are low. This blended sequencing technique also offers a low cost protocol in developing countries like Bangladesh.

Genomic variability correlates with biofilm phenotypes in multidrug resistant clinical isolates of Pseudomonas aeruginosa.

The multifactorial nature of Pseudomonas aeruginosa biofilm development and genomic variabilities implicates its resistance to conventional antimicrobials and virulence. Therefore, genetic determinants need to be extensively studied to block the early steps of biofilm or already formed biofilms. In this study, a total of 20 multidrug resistant (MDR) clinical P. aeruginosa isolates were evaluated for their biofilm forming abilities and related genes. Of the isolates tested, all of them showed surface attachment tendencies in nutrient limiting conditions, and classified as strong (SBF = 45%), moderate (MBF = 30%) and weak (WBF = 25%) biofilm formers. Complete genome sequencing of representative strong (DMC-27b), moderate (DMC-20c) and weak biofilm former (DMC-30b) isolates was performed. Analysis of biofilm related genes in the sequenced genomes revealed that, 80 of the 88 biofilm related genes possess 98-100% sequence identity to the reference PAO1 strain. Complete and partial sequence data of LecB proteins from tested isolates indicate that isolates containing PA14-like LecB sequences produced strong biofilms. All of the 7 pel operon protein coding genes in weak biofilm former isolate 30b showed significant nucleotide sequence variation with other tested isolates, and their corresponding proteins are 99% identical with the pel operon proteins of PA7. Bioinformatics analyses identified divergent sequence and structural features that separate PA7 like pel operon proteins from reference PAO1-like pel operon. Congo red and pellicle forming assays revealed that the sequence and structure variations may have interfered with the Pel production pathway and resulted in impaired Pel production in isolate 30b that has a PA7 like pel operon. Expression analysis also showed that both pelB and lecB genes were about 5 to 6 folds upregulated after 24 h in SBF 27b in comparison with WBF 30b. Our findings indicate significant genomic divergence in biofilm related genes of P. aeruginosa strains that affect their biofilm phenotypes.

Temporal dynamics and fatality of SARS-CoV-2 variants in Bangladesh.

Background and Aims: Since the beginning of the SARS-CoV-2 pandemic, multiple new variants have emerged posing an increased risk to global public health. This study aimed to investigate SARS-CoV-2 variants, their temporal dynamics, infection rate (IFR) and case fatality rate (CFR) in Bangladesh by analyzing the published genomes. Methods: We retrieved 6610 complete whole genome sequences of the SARS-CoV-2 from the GISAID (Global Initiative on Sharing all Influenza Data) platform from March 2020 to October 2022, and performed different in-silico bioinformatics analyses. The clade and Pango lineages were assigned by using Nextclade v2.8.1. SARS-CoV-2 infections and fatality data were collected from the Institute of Epidemiology Disease Control and Research (IEDCR), Bangladesh. The average IFR was calculated from the monthly COVID-19 cases and population size while average CFR was calculated from the number of monthly deaths and number of confirmed COVID-19 cases. Results: SARS-CoV-2 first emerged in Bangladesh on March 3, 2020 and created three pandemic waves so far. The phylogenetic analysis revealed multiple introductions of SARS-CoV-2 variant(s) into Bangladesh with at least 22 Nextstrain clades and 107 Pangolin lineages with respect to the SARS-CoV-2 reference genome of Wuhan/Hu-1/2019. The Delta variant was detected as the most predominant (48.06%) variant followed by Omicron (27.88%), Beta (7.65%), Alpha (1.56%), Eta (0.33%) and Gamma (0.03%) variant. The overall IFR and CFR from circulating variants were 13.59% and 1.45%, respectively. A time-dependent monthly analysis showed significant variations in the IFR (p = 0.012, Kruskal-Wallis test) and CFR (p = 0.032, Kruskal-Wallis test) throughout the study period. We found the highest IFR (14.35%) in 2020 while Delta (20A) and Beta (20H) variants were circulating in Bangladesh. Remarkably, the highest CFR (1.91%) from SARS-CoV-2 variants was recorded in 2021. Conclusion: Our findings highlight the importance of genomic surveillance for careful monitoring of variants of concern emergence to interpret correctly their relative IFR and CFR, and thus, for implementation of strengthened public health and social measures to control the spread of the virus. Furthermore, the results of the present study may provide important context for sequence-based inference in SARS-CoV-2 variant(s) evolution and clinical epidemiology beyond Bangladesh.

Contamination of Selected Toxic Elements in Integrated Chicken-Fish Farm Settings of Bangladesh and Associated Human Health Risk Assessments.

The presence of trace elements in the environment can contaminate a food chain of an agro farm in various ways. Integrated chicken-fish farms (i.e., where poultry chicken and fish are cultivated in same places) are getting popular nowadays to meet the demands of a balanced diet. The present study conducted a health risk assessment on the basis of selected heavy metal (i.e., Cr and Pb) and metalloid (i.e., As) contamination in this type of farm in Bangladesh. Samples of various types were collected from different farms between September 2019 and March 2020. The concentrations of the elements were checked by Flame-AAS and HG-AAS. Our findings demonstrated that the elements' concentrations in fishes were simultaneously induced by the habitation and bioaccumulation through the food chain of the farm. The concentrations of As and Pb in the chicken parts and Cr and As in some fishes were greater than the highest limits set by different permissible standards. Overall, the metal concentration obtained in different samples was in descending order: sediment > droppings > different fish parts > various species of chicken > pond water. Among the pollutants, As gave target hazard quotient (THQ) values higher than 1 for all the species, suggesting health risks from the intake of fishes and chicken. However, there was non-target cancer risk present while considering all the elements together. Notably, the study found carcinogenic risks of As, Pb, and Cr for humans due to poultry and/or fish consumption; the identified health risks associated with the integrated farming setting will be crucial in further tackling strategies. Investigation of the possible sources of heavy metals in commercial chicken feeds and regular monitoring of groundwater used for agro-farming are highly recommended to reduce the burden.

Microbiome analysis revealing microbial interactions and secondary bacterial infections in COVID-19 patients comorbidly affected by Type 2 diabetes.

The mortality of coronavirus disease 2019 (COVID-19) disease is very high among the elderly or individuals having comorbidities such as obesity, cardiovascular diseases, lung infections, hypertension, and/or diabetes. Our study characterizes the metagenomic features in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-infected patients with or without type 2 diabetes, to identify the microbial interactions associated with its fatal consequences.This study compared the baseline nasopharyngeal microbiome of SARS-CoV-2-infected diabetic and nondiabetic patients with controls adjusted for age and gender. The metagenomics based on next-generation sequencing was performed using Ion GeneStudio S5 Series and the data were analyzed by the Vegan-package in R. All three groups possessed significant bacterial diversity and dissimilarity indexes (p < 0.05). Spearman's correlation coefficient network analysis illustrated 183 significant positive correlations and 13 negative correlations of pathogenic bacteria (r = 0.6-1.0, p < 0.05), and 109 positive correlations between normal flora and probiotic bacteria (r > 0.6, p < 0.05). The SARS-CoV-2 diabetic group exhibited a significant increase in pathogens and secondary infection-causing bacteria (p < 0.05) with a simultaneous decrease of normal flora (p < 0.05). The dysbiosis of the bacterial community might be linked with severe consequences of COVID-19-infected diabetic patients, although a few probiotic strains inhibited numerous pathogens in the same pathological niches. This study suggested that the promotion of normal flora and probiotics through dietary supplementation and excessive inflammation reduction by preventing secondary infections might lead to a better outcome for those comorbid patients.

Genomic analysis and in vivo efficacy of Pediococcus acidilactici as a potential probiotic to prevent hyperglycemia, hypercholesterolemia and gastrointestinal infections.

Lactic acid bacteria are the well acknowledged probiotics that can cure a variety of diseases. In this study, we observed the in vivo potentials of Pediococcus to treat hyperglycemia, hypercholesterolemia and gastrointestinal infections. A total of 77 Lactobacillus were isolated from the milk of 10 cows and 10 goats, four of those strains inhibited both carbohydrates-hydrolyzing enzymes, α-glucosidase, and α-amylase. They all showed antagonistic effects on pathogenic E. coli and S. Typhimurium which were confirmed by performing pathogen challenge test and visualizing on Electron microscopy. 16S rRNA gene sequence identified that all four strains belong to Pediococcus genus which were further distinguished as Pediococcus acidilactici by pheS gene sequence. Whole genome sequence analysis revealed their non-pathogenic properties for human and the presence of probiotic genes responsible for stress resistance, immunomodulation, adhesion, metal and drug resistance. In vivo trial with diabetes-induced mice ascertained that all Pediococcus acidilactici had significant potentials to reduce elevated glucose and low-density lipoprotein level in blood. Interestingly, two out of four strains were significantly more effective (p < 0.0001 each) than metformin in reducing the blood glucose level. This in vivo study demonstrated that Pediococcus acidilactici might be a promising probiotic to prevent hyperglycemia, hypercholesterolemia and gastrointestinal infections.

Genomic analysis of SARS-CoV-2 variants of concern identified from the ChAdOx1 nCoV-19 immunized patients from Southwest part of Bangladesh.

BACKGROUND: Bangladesh introduced ChAdOx1 nCoV-19 since February, 2021 and in six months, only a small population (12.8%) received either one or two dose of vaccination like other low-income countries. The COVID-19 infections were continued to roll all over the places although the information on genomic variations of SARS-CoV-2 between both immunized and unimmunized group was unavailable. The objective of this study was to compare the proportion of immune escaping variants between those groups. METHODS: A total of 4718 nasopharygeal samples were collected from March 1 until April 15, 2021, of which, 834 (18%) were SARS-CoV-2 positive. The minimum sample size was calculated as 108 who were randomly selected for telephone interview and provided consent. The prevalence of SARS-CoV-2 variants and disease severity among both immunized and unimmunized groups was measured. A total of 63 spike protein sequences and 14 whole-genome sequences were performed from both groups and phylogenetic reconstruction and mutation analysis were compared. RESULTS: A total of 40 respondents (37%, N = 108) received single-dose and 2 (2%) received both doses of ChAdOx1 nCoV-19 vaccine, which significantly reduce dry cough, loss of appetite and difficulties in breathing compared to none. There was no significant difference in hospitalization, duration of hospitalization or reduction of other symptoms like running nose, muscle pain, shortness of breathing or generalized weakness between immunized and unimmunized groups. Spike protein sequence assumed 21 (87.5%) B.1.351, one B.1.526 and two 20B variants in immunized group compared to 27 (69%) B.1.351, 5 (13%) B.1.1.7, 4 (10%) 20B, 2 B.1.526 and one B.1.427 variant in unimmunized group. Whole genome sequence analysis of 14 cases identified seven B.1.351 Beta V2, three B.1.1.7 Alpha V1, one B.1.526 Eta and the rest three 20B variants. CONCLUSION: Our study observed that ChAdOx1 could not prevent the new infection or severe COVID-19 disease outcome with single dose while the infections were mostly caused by B.1.351 variants in Bangladesh.

Dominant clade-featured SARS-CoV-2 co-occurring mutations reveal plausible epistasis: An in silico based hypothetical model.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved into eight fundamental clades with four of these clades (G, GH, GR, and GV) globally prevalent in 2020. To explain plausible epistatic effects of the signature co-occurring mutations of these circulating clades on viral replication and transmission fitness, we proposed a hypothetical model using in silico approach. Molecular docking and dynamics analyses showed the higher infectiousness of a spike mutant through more favorable binding of G614 with the elastase-2. RdRp mutation p.P323L significantly increased genome-wide mutations (p < 0.0001), allowing for more flexible RdRp (mutated)-NSP8 interaction that may accelerate replication. Superior RNA stability and structural variation at NSP3:C241T might impact protein, RNA interactions, or both. Another silent 5'-UTR:C241T mutation might affect translational efficiency and viral packaging. These four G-clade-featured co-occurring mutations might increase viral replication. Sentinel GH-clade ORF3a:p.Q57H variants constricted the ion-channel through intertransmembrane-domain interaction of cysteine(C81)-histidine(H57). The GR-clade N:p.RG203-204KR would stabilize RNA interaction by a more flexible and hypo-phosphorylated SR-rich region. GV-clade viruses seemingly gained the evolutionary advantage of the confounding factors; nevertheless, N:p.A220V might modulate RNA binding with no phenotypic effect. Our hypothetical model needs further retrospective and prospective studies to understand detailed molecular events and their relationship to the fitness of SARS-CoV-2.

A rapid and cost-effective multiplex ARMS-PCR method for the simultaneous genotyping of the circulating SARS-CoV-2 phylogenetic clades.

Tracing the globally circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) phylogenetic clades by high-throughput sequencing is costly, time-consuming, and labor-intensive. We here propose a rapid, simple, and cost-effective amplification refractory mutation system (ARMS)-based multiplex reverse-transcription polymerase chain reaction (PCR) assay to identify six distinct phylogenetic clades: S, L, V, G, GH, and GR. Our multiplex PCR is designed in a mutually exclusive way to identify V-S and G-GH-GR clade variants separately. The pentaplex assay included all five variants and the quadruplex comprised of the triplex variants alongside either V or S clade mutations that created two separate subsets. The procedure was optimized with 0.2-0.6 µM primer concentration, 56-60°C annealing temperature, and 3-5 ng/µl complementary DNA to validate on 24 COVID-19-positive samples. Targeted Sanger sequencing further confirmed the presence of the clade-featured mutations with another set of primers. This multiplex ARMS-PCR assay is a fast, low-cost alternative and convenient to discriminate the circulating phylogenetic clades of SARS-CoV-2.

Emergence of European and North American mutant variants of SARS-CoV-2 in South-East Asia.

The SARS-CoV-2 coronavirus is responsible for the current COVID-19 pandemic, with an ongoing toll of over 5 million infections and 333 thousand deaths worldwide within the first 5 months. Insight into the phylodynamics and mutation variants of this virus is vital to understanding the nature of its spread in different climate conditions. The incidence rate of COVID-19 is increasing at an alarming pace within subtropical South-East Asian nations with high temperatures and humidity. To understand this spread, we analysed 444 genome sequences of SARS-CoV-2 available on the GISAID platform from six South-East Asian countries. Multiple sequence alignments and maximum-likelihood phylogenetic analyses were performed to analyse and characterize the non-synonymous (NS) mutant variants circulating in this region. Global mutation distribution analysis showed that the majority of the mutations found in this region are also prevalent in Europe and North America, and the concurrent presence of these mutations at a high frequency in other countries indicates possible transmission routes. Unique spike protein and non-structural protein mutations were observed circulating within confined area of a given country. We divided the circulating viral strains into four major groups and three subgroups on the basis of the most frequent NS mutations. Strains with a unique set of four co-evolving mutations were found to be circulating at a high frequency within India, specifically. Group 2 strains characterized by two co-evolving NS mutants which alter in RdRp (P323L) and spike (S) protein (D614G) were found to be common in Europe and North America. These European and North American variants have rapidly emerged as dominant strains within South-East Asia, increasing from a 0% prevalence in January to an 81% by May 2020. These variants may have an evolutionary advantage over their ancestral types and could present a large threat to South-East Asia for the coming winter.

Understanding the possible origin and genotyping of the first Bangladeshi SARS-CoV-2 strain.

The novel coronavirus, SARS-CoV-2, has caused the most unfathomable pandemic in the history of humankind. Bangladesh is also a victim of this critical situation. To investigate the genomic features of the pathogen from Bangladesh, the first complete genome of the virus has very recently been published. Therefore, long-awaited questions regarding the possible origin and typing of the strain(s) can now be answered. Here, we endeavor to mainly discuss the published reports or online-accessed data (results) regarding those issues and present a comprehensive picture of the typing of the virus alongside the probable origin of the subclade containing the Bangladeshi strain. Our observation suggested that this strain might have originated from the United Kingdom or the other European countries epidemiologically linked to the United Kingdom. According to different genotyping classification schemes, this strain belongs to the A2a clade under the G major clade, is of B and/or L type, and is a SARS-CoV-2a substrain. In the future, randomized genomic data will certainly increase in Bangladesh, however because of globalization and immigrant movement, we urgently need a mass regional sequencing approach targeting the partial or complete genome that can link the epidemiological data and may help in further clinical intervention.

Multidrug-Resistant Avian Pathogenic Escherichia coli Strains and Association of Their Virulence Genes in Bangladesh.

The avian pathogenic Escherichia coli (APEC) strains are the chief etiology of colibacillosis worldwide. The present study investigated the circulating phylotypes, existence of virulence genes (VGs), and antimicrobial resistance (AMR) in 392 APEC isolates, obtained from 130 samples belonged to six farms using both phenotypic and PCR-based molecular approaches. Congo red binding (CRB) assay confirmed 174 APEC isolates which were segregated into ten, nine, and eight distinct genotypes by RAPD assay (discriminatory index, DI = 0.8707), BOX-PCR (DI = 0.8591) and ERIC-PCR (DI = 0.8371), respectively. The combination of three phylogenetic markers (chuA, yjaA and DNA fragment TspE4.C2) classified APEC isolates into B23 (37.36%), A1 (33.91%), D2 (11.49%), B22 (9.20%), and B1 (8.05%) phylotypes. Majority of the APEC isolates (75-100%) harbored VGs (ial, fimH, crl, papC, and cjrC). These VGs (papC and cjrC) and phylotypes (D2 and B2) of APEC had significant (p = 0.004) association with colibacillosis. Phylogenetic analysis showed two distinct clades (clade A and clade B) of APEC, where clade A had 98-100% similarity with E. coli APEC O78 and E. coli EHEC strains, and clade B had closest relationship with E. coli O169:H41 strain. Interestingly, phylogroups B2 and D2 were found in the APEC strains of both clades, while the strains from phylogroups A1 and B1 were found in clade A only. In this study, 81.71% of the isolates were biofilm formers, and possessed plasmids of varying ranges (1.0 to 54 kb). In vitro antibiogram profiling revealed that 100% isolates were resistant to ≥3 antibiotics, of which 61.96%, 55.24%, 53.85%, 51.16% and 45.58% isolates in phylotypes B1, D2, B22, B23, and A1, respectively, were resistant to these antimicrobials. The resistance patterns varied among different phylotypes, notably in phylotype B22, showing the highest resistance to ampicillin (90.91%), nalidixic acid (90.11%), tetracycline (83.72%), and nitrofurantoin (65.12%). Correspondence analysis also showed significant correlation among phylotypes with CRB (p = 0.008), biofilm formation (p = 0.02), drug resistance (p = 0.03), and VGs (p = 0.06). This report demonstrated that B2 and A1 phylotypes are dominantly circulating APEC phylotypes in Bangladesh; however, B2 and D2 are strongly associated with the pathogenicity. A high prevalence of antibiotic-resistant APEC strains from different phylotypes suggest the use of organic antimicrobial compounds, and/or metals, and the rotational use of antibiotics in poultry farms in Bangladesh.

Insights Into the Resistome of Bovine Clinical Mastitis Microbiome, a Key Factor in Disease Complication.

Bovine clinical mastitis (CM) is one of the most prevalent diseases caused by a wide range of resident microbes. The emergence of antimicrobial resistance in CM bacteria is well-known, however, the genomic resistance composition (the resistome) at the microbiome-level is not well characterized. In this study, we applied whole metagenome sequencing (WMS) to characterize the resistome of the CM microbiome, focusing on antibiotics and metals resistance, biofilm formation (BF), and quorum sensing (QS) along with in vitro resistance assays of six selected pathogens isolated from the same CM samples. The WMS generated an average of 21.13 million reads (post-processing) from 25 CM samples that mapped to 519 bacterial strains, of which 30.06% were previously unreported. We found a significant (P = 0.001) association between the resistomes and microbiome composition with no association with cattle breed, despite significant differences in microbiome diversity among breeds. The in vitro investigation determined that 76.2% of six selected pathogens considered "biofilm formers" actually formed biofilms and were also highly resistant to tetracycline, doxycycline, nalidixic acid, ampicillin, and chloramphenicol and remained sensitive to metals (Cr, Co, Ni, Cu, Zn) at varying concentrations. We also found bacterial flagellar movement and chemotaxis, regulation and cell signaling, and oxidative stress to be significantly associated with the pathophysiology of CM. Thus, identifying CM microbiomes, and analyzing their resistomes and genomic potentials will help improve the optimization of therapeutic schemes involving antibiotics and/or metals usage in the prevention and control of bovine CM.